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1 Abstract 



■ The Crab Nebula is well-described at optical wavelengths and above by a 
steady-state magnetohydrodynamic model due to Kennel and Coroniti (1984). 
Can this class of model describe other pulsar-wind nebulae? I exhibit simple 
generalizations of KC models, for various values of a, the ratio of magnetic to 
particle flux input at the wind shock. I calculate the evolution of the electron 

^1 , spectrum and synchrotron emissivity in the nebula assuming spherical sym- 

00 ' metry, a steady state, and a purely toroidal magnetic field. Emission profiles 

, depend on the initial magnetic field Bq, the electron spectral index s, and the 

■ angle of the toroidal axis with the line of sight, <p. I show integrated spectra 
and radial profiles for various cases, along with predicted variations of pho- 

' ton index F with radius in X-rays. Most models predict much smaller sizes 

in X-rays, and curves of r{r) which do not resemble observations. Further 
Q ■ elaborations of the dynamics of PWNe seem necessary. 

! 2 Introduction 



Pulsar- wind nebulae (PWNe) are bubbles of relativistic particles and magnetic 
5^ ■ field inflated by a pulsar, forming center-brightened radio and X-ray sources 

I of synchrotron radiation , [2] , [HI , 0] . Most early work on PWNe assumed 

they were surrounded by shell supernova remnants (SNRs), seen or unseen, 
and interacted with the interior of those remnants. Before the passage of a 
reverse shock, the PWN would expand into unshocked ejecta in the SNR 
interior; afterwards, it would be compressed by the reverse-shock passage and 
expand much more slowly into the thermalized, shocked ejecta [2]. 

The most elaborate early model for the high-energy (optical and above) 
emission from PWNe was that of Kennel & Coroniti [3], 0] who solved the 
steady-state relativistic MHD equations in spherical symmetry with an az- 
imuthal magnetic field to obtain velocity and density profiles of the outflow- 
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ing relativistic material. They showed that a strong (high compression) MHD 
shock required a fluid dominated by particles, not magnetic fleld, as parame- 
terized by cr = {B / An) / n'j'^ Pmc^ , the ratio of magnetic (plus electric) energy 
flux to particle flux just upstream of the wind shock. Here n is the comoving 
density and 7 the flow Lorentz factor. Applying the Rankine-Hugoniot jump 
conditions at the wind shock (radius r^) showed that unless cr <C 1, the post- 
shock fluid would remain at high speed, within a factor of 3 of the upstream 
(relativistic) speed - clearly inconsistent with the observed expansion of the 
Crab Nebula at w ~ 2000 km s~^. Behind the shock, pure flux- freezing meant 
that the toroidal magnetic field strength evolved as B oc pr, with p the mat- 
ter density. Kennel & Coroniti ^ were able then to follow the evolution of 
electron energies and predict both the integrated spectrum and the brightness 
profile of remnants. They specialized to the Crab Nebula and showed calcula- 
tions primarily for that case. They were able to describe the optical through 
hard X-ray spectrum well, as well as the gross variation of nebular size with 
X-ray energy ("size" being typically the FWHM of emission measured with 
relatively crude imaging instruments). However, the high value of preshock 
Lorentz factor of the wind they found (7 ^ 10^) made it difficult to explain 
the radio emission which requires much lower energies, and more electrons, 
than would be produced by thermalizing this extremely relativistic wind. 

Since the advent of the newer generation of X-ray telescopes, beginning 
with ASCA and continuing with the Chandra X-ray Observatory and XMM- 
Newton, we have seen an explosion in the number and variety of PWNe in 
X-rays, and have observed a much wider range of morphological type. It is 
of interest to ask if the class of steady-state MHD models proposed by KC 
can be generalized appropriately to describe other PWNe, and if not, to get 
an idea of what kind of extensions might be necessary. This is the project 
described in this paper. 

One major problem in modeling several larger, fainter PWNe such as 3C 
58 |H] appears in comparing radio and X-ray images. In Kennel & Coroniti's 
model for the Crab, synchrotron losses on outflowing particles cause the neb- 
ular size to decrease with frequency at optical and higher frequencies. In fact 
the Crab is considerably smaller at optical than radio frequencies, and contin- 
ues to diminish in size roughly as i? oc 1/^0.15 j^gy However, 3C 
58 shows X-ray emission extending to near the edges of the radio emission, 
though it is faint there (Slane, private communication). These results call into 
question the MHD picture of transport of particle energy and magnetic flux 
throughout the volume of a PWN. 

It is straightforward to show from the MHD equations |2j that postshock 
flow solutions have accelerating and decelerating branches, with the latter 
beginning at vq = c/3 and initially decelerating as (pure hydrodynamic, 
isobaric flow; magnetic energy is unimportant since a must be small to fit the 
outer nebular boundary conditions). However, in this case B cc pr oa 1/rv r, 
so the magnetic energy rises until it becomes dynamically important, slowing 
the deceleration and causing the velocity to approach a constant value, with 



Magnetohydrodynamic Models for the Structure of Pulsar- Wind Nebulae 3 

B l/r. The full hydrodynamic solution can be well approximated by a 
simple analytic form: defining u = v/vq, 

U{r) = ^{1 - Uoo) + Uoo, (1) 

where the asymptotic velocity is Voc/c = cr/(l + <j) (Fig. 1). 

Approximate PWN Velocity Profiles Approximate PWN Magnetic-Field Profiles 
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Fig. 1. Solid lines: a = 0.001. Dashed lines: a = 0.0032. Dotted lines: a = 0.01. 
Dot-dashed lines: a — 0.032. 



3 Particle evolution and luminosity calculation 

Electron energies E evolve due both to radiative and adiabatic losses: E = 
{dE/dV){dV/dt) — aB'^E^, where for synchrotron losses a = 1.57 x 10^'^ cgs. 
This equation can be integrated to give j^, [7] 

^ Epa^/^ ^ Eoa^ 

1 + aEojB^a^/^dt i + Eo J z-^u'^^/^dz 

where a = p/po, z = r/rp, and u = v/vq. The quantity Ef = {aB^r^ / vq)~^ is 
a "fiducial energy" , that an initially infinitely energetic electron would have 
after radiating in a field Bq for a time tq/vq (basically a parameterization 
of -Bo)- From this result, the evolution of an arbitrary distribution can be 
calculated: 

. mE,/JI'^ . NiE,iE„'^. . NiEm> (I) {^)"' 

Given the approximate dynamics above, the electron distribution and syn- 
chrotron emissivity can be calculated numerically at each radius. If magnetic 
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Electron Spectrum Inside PWN 



Synchrotron Emissivity Inside PWN 
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Fig. 2. Flat-spGctrum, high-ir case: s = 1.0, a = 0.01, Ef = 910 erg, and Emin = 
10~^ erg. From left to right, lines show the spectrum and emissivity at five radii 
(outside in): r/vo = 10, 8, 6, 4, 2. 



PWN Synchrotron Spectrum 



Radial Variation of PWN Photon Index (1-10 keV) 
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Dashed hne: Intensity Profile at 1 keV 
Dot— dash line: Intensity Profile at 10 keV 



Fig. 3. Same parameters as above. Note curvature of the integrated spectrum due 
to the energy- loss "bump" . The solid line on the right is the 1-10 keV photon index; 
the profiles are for aspect angle = 0° (B in the plane of the sky). 



PWN Synchrotron Spectrum 



Radial Variation of PWN Photon Index (1-10 keV) 
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Dashed line: Intensity Profile at 1 keV 
Dot— dash line; Intensity Profile at 10 keV 



Fig. 4. Steeper-spectrum, low-a case: s = 1.5, a = 0.001, Ef = 9.1 x 10^ erg, and 

aspect angle = 90°. 
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field lines are circles in an equatorial plane, an elementary rotation can give 
brightness profiles as a function of aspect angle (p between the polar axis and 
the line of sight. The synchrotron emissivity was integrated numerically to 
generate the brightness profiles shown above, for an injected flat power-law 
N{E) = KE~^ with s < 2 as deduced from observed PWN radio spectra. 
The flat spectra cause a "bump" to appear just below the high-energy cutoff 
that appears as a result of radiative losses. 

4 Conclusions 

Steady-state MHD models of pulsar-wind nebulae can be well approximated 
by a simple velocity law, which then dictates the evolution of the injected 
electron distribution and the integrated spectrum and proflle of emission, 
depending on a. Spectral properties also depend on B and spectral index s. 
Varying the aspect angle changes the flux normalization but has no effect on 
the integrated spectrum and little on the proflle shape at photon energies for 
which losses are important. For very flat injected spectra (s ~ 1), the effects 
of the energy-loss "bump" are signiflcant, causing concave- up curvature over 
a broad frequency range. 

The models are in signiflcant conflict with the observed radial dependence 
of photon index F in X-rays (e.g., G21.5-0.9 8J, which appears roughly linear 
with r. All models drop in size by at least a factor 2 between radio and X-rays, 
but some observed PWNe do not show this. Unless magnetic-field strengths 
are unrealistically low, most X-ray profiles have HWHM radii only 3-6 times 
the injection radius - also perhaps at odds with observations. 

We conclude that there are significant discrepancies between the predic- 
tions of the simple MHD models and X-ray observations of PWNe. Nonsteady 
or nonspherical flows, and/or electron transport by diffusion as well as con- 
vection (e.g., i^,), may be necessary to explain the observations. 
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